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Abstract. Wc investigate spin dynamics of a two-component Bose-Einstein 
condensates with weak Josephson coupHng. Analytical expressions of atom- 
number squeezing and bipartite entanglement are presented for atom-atom repulsive 
interactions. For attractive interactions, there is no number squeezing; however, the 
squeezing parameter is still useful to recognize the appearance of Schrodinger's cat 
state. 
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Atom-number squeezing has attracted much attention due to its potential 
apphcations in quantum metrology and quantum information [B [21 [3l HI O [6]. As 
a special case of spin squeezed states [71 El [9], [101 El [El [iSl El [15], the number 
squeezed state shows reduced spin fluctuation of the component below the standard 
quantum limit (SQL), which in turn suppresses the deleterious effects arising from phase 
diffusion [16l [T71 [HI [191 [20l [21] . Dynamical generation of the number squeezing has been 
theoretically investigated based upon a two-mode boson model with a weak Josephson 
couphng [221 [231 [251 [261 [271 [28] . 

So far, the number squeezing has been demonstrated indirectly by detecting an 
increased phase fluctuation [Il[6], but not atom-number variance. However, the variances 
of relative number and phase has a nontrivial relation so a direct measurement of 
the number fluctuation is of interests and necessary [27]. In this letter, we study 
atom-number squeezing and bipartite entanglement of a two-component BEC with the 
Josephson coupling. Two analytical expressions are discovered, which provides us a 
more direct way to measure the number variance and entropy of entanglement through 
extracting phase coherence (i.e., the visibility) in atomic interference experiments. 

We consider a two-component Bose-Einstein condensate with hyperfine states |1) 
and 1 2) coupled by an external microwave (or rf) field [291 [201 [SH [32]. For the 
BEC confined in a deep three-dimensional harmonic potential, we adopt single-mode 
approximation [331 [311 [Ml [36], which results in the following Hamiltonian {h = 1): 

H = 6J,- n.h + xJl (1) 

where the detuning 5, the Josephson-like coupling constant fi, and the mean-field 
interaction strength x can be controlled artificially. Angular momentum operators 
J+ = {J-Y = a\ai and = {N2 — Ni)/2 are introduced, where and N„ (= aJ.ao-) 
are the annihilation and number operators for the two modes a = 1, 2. Eq. ([T]) is 
equivalent with a two-site Bose-Hubbard hamiltonian [37], where 6 denotes potential 
bias of a double well, Q the hopping term, and x the on-site interaction. Total particle 
number = A'^i -|- A^2 is assumed to be a fixed c number, so atomic number operators 
= j + {—lyjz with 0" = 1,2 and j = N/2. Atom number fluctuations, defined as 
usual {AN^y = {N^) — {N„)^ are the same and equal to the variance (AJ^)^. 

For any state vector |\E'), one can determine the mean spin: (J) = {{Jx), {Jy)i {Jz)), 
where {Jx) = Re(J+) and {Jy) = Im(J+). The expectation value (J+) relates to the 
first-order cross correlation function [201 [21] : 

(1) ^ I (4^1 ) I ^ /2^ 

which measures phase coherence of the two-component BEC. It is observable by 
extracting the visibility of atomic interference fringes [H [21 [3l [H [5l [6l [311 [371 [21] . Similar 
definition of the phase coherence has been proposed in Refs. [39l [ID]. The degree of 
atom-number squeezing is quantified by a parameter [H [22] 

2 ^ A^(AA^.)^ ^ 2j(AJ.)^ . . 

^ {N,){N2) f-{J.r 
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where total particle number = 2j is assumed to be fixed. It has been shown that 
there is no squeezing (i.e., = 1) for coherent spin state (CSS) [H]: 

\0,(f)) = exp[i9 {Jr, sin (j) - Jy cos (f))]\j , j) , (4) 

which yields (J^) = jcos9 and (AJ^)^ = (j'/2)sin^^. The condition of the number 
squeezing is therefore < 1 [22], which is consistent with previous one [6]: (A J^)^ < j/2 
for a symmetric BEC with population imbalance (Jz) = 0. 

Firstly, let us consider spin dynamic of the symmetric BEC {5 = 0) with repulsive 
interactions (x > 0) for an initial CSS state |\E'(0)) = |f,0) = \j,j)x, which is also an 
eigenvector of with eigenvalue j. Experimentally, such a state has been prepared by 
applying a two-photon 7i/2 pulse to the condensed atoms occupied in the internal state 
|2) [29]. The state vector at any time t can be expanded as |\E') = J2m^rn\j,'m'), where 
the probability amplitudes are determined by time-dependent Schrodinger equation 
with the initial condition: Cm(0) = {j,'m\j,j)x = ^(j+m)^^^- Note that the initial state 
shows the population imbalance (Jz(0)) = and the number variance {AJz{0)y = j/2 

(i.e., e = 1). 

As the simplest case, we consider the Hamiltonian ([1]) with = 0. This is the 
one-axis twisting model, proposed originally by Kitagawa and Ueda [7]. The phase 
coherence can be solved exactly as g^^\t) = cos^''~^(x^)- In the short-time limit, it decays 
exponentially as g^^\t) ~ exp[— (t/t^)^] with x^d = denoting a characteristic time 

scale for phase coherence [Ml [171 dH [19] . The damping of phase coherence, known as 
phase diffusion [161 (El O UHl 1201 [21] has been observed in experiment by extracting 
the visibility of the Ramsey fringe [38]. During phase diffusion, the number squeezing 
remains constant because of conserved probability distribution |cmP- 

The interplay between the nonlinear interaction and the Josephson coupling leads 
to suppressed phase diffusion due to the appearance of number squeezing. So far, the 
number squeezing has been demonstrated for the BEC in optical lattices [B 121 S H], 
optical trap [5], and atom chip [6]. To understand how this works, we consider 
unitary evolution of the symmetric BEC governed by Hamiltonian ([T|). In Schrodinger 
picture, the amplitudes always satisfy the relation c_m = Cm, which in turn leads to 
{Jy) = {Jz) = and {Jx) = {J+) 7^ 0. Moreover, Eq. ([2]) and Eq. ([3]) now reduce to 
g^^^ = \{Jx)\/j and = 2(J^)/j. A relation between ^f^-^^ and ^'^ can be obtained by 
examining Heisenberg equations of motion: 

Jx = - X{JyJz + JzJy), (5) 
jy = QJz + X{JxJz + JzJx), (6) 
jz= - ^Jy. (7) 

For time-independent Q and x, we have dJ^/dt = {Q/x)dJx/dt, and thus 

e{t) = 1-^-^ [IT (8) 

with the upper sign for the mean spin {Jx) > and the lower one for {Jx) < 0. 
Eq. ([8]) provides us an exact relation between the number squeezing and the 
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Figure 1. (Color online) Time evolution of phase coherence g*^^^ (black thin) and 
number-squeezing parameter (red) for N = 20, and the coupling: ri/|x| = 1.732 (a) 
and (c), r^/lxl = 10 (b) and (d). In (c) and (d), negative x case is considered. The 
circles and the dashed green line are given by Eq. ([5]) and Eq. respectively. The 
time is in units of 

phase coherence g^^\ To confirm it, we consider two-particle (A^ = 2) case. It 
is easy to obtain g^^'^ = 1 — |(xsina;2V^2)^ and = 1 — r2x(sinu;2V'^2)^) with 
LU2 = (fi^ + Obviously, Eq.(l8]) holds for the two-particle case. In addition, 

we find that local minimum of ^"^ (i.e., maximal number squeezing) occurs at time 
^min = 7i"/(2a;2). If an optimal coupling Q/x = 1/2 is applied, the system will evolve 
into maximally number-squeezed state (MSS): |\1/)mss = 'ie~*'^/(^^''|l, 0), which exhibits 
perfect squeezing (AJj)^ = [8]. 

For large case, there exist no exact solutions; however, some approximated 
solutions are obtainable if the coupling is strong enough. In this case, the mean spin 
(Jx) almost remains unchanged at j. As a result, we adopt frozen-spin approximation 
(FSA) P, [221 El], i-e., replacing by j (= A^/2) in the Heisenberg equations, and obtain 
Jz = —Vl^Jz — ^x{JxJz + JzJx) — —oo%Jz, where 00^ = >y^{^ + A^x) [23]. Solving the 
above equation, it is easy to get harmonic solutions of Jy{t) and Jz{t) [21]. Inserting 
them into Eq. ([5]), we obtain Ja.(t) and also (7*^^^ = \ {Jx)\/j — 1 — Y(xsinu;Art/a;Ar)^. Now, 
the phase coherence g^^\t) almost exhibits sinusoidal oscillations [thin lines of Fig. [1], 
but not exponential damping as previous case. In other words, the phase diffusion is 
strongly suppressed due to the Josephson coupling. From Eq. (|8]), we further obtain 

em.l-N(lx(?^)\ (9) 

\ UJn J 

Clearly, the maximal squeezing appears at time tmm = tt/ {2u!n) with, ^^j^^ = fl'^/uj^ < 1 
[231 [2l] . Dunningham et al. [25] has independently derived the time tmin using a 
semiclassical analysis of Hamiltonian ([1]) [3^. For large A^ case (> 10'^), the optimal 
coupling obeys the power rule Q/x ~ 0.58A^°'^^ [261 [2Z]- In Fig. i^a), time evolution 
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of 5f(^) and is plotted for the optimal coupling Q/x = 1.732 and = 20. The FSA 
works (green dashed curves) quite well to predict tmin. For relatively large fl/X: say 
n/x = 10 for N = 20, the FSA follows full evolution of [see Fig. W[h)]. As shown 
Fig. [T](c) and (d), Eq.® keeps hold for negative x case [see below]. 

The relation between spin squeezing and quantum entanglement is of interests. 
It was shown that the obtained squeezing is useful for quantum metrology [8], [9] and 
entanglement [12], provided that Cs = V^i^Jz)/\{J+)\ = i/g''^^ < 1 [1]. For the system 
considered here, however, only entanglement between the two modes (i.e., bipartite 
entanglement) is accessible due to the indistinguishability of identical bosons [43]. A 
standard measure of bipartite entanglement is the so-called entropy of entanglement 

[121 m]: 

E{t) = - Tr[pi log(pi)] = -Y, log (|c^n , (10) 

m,=-j 

where pi = Tr2(p) is the reduced density operator for mode 1 obtained by partial trace 
over mode 2. The value of E varies between 0, for the separable product states to a 
maximum of \og{d), for maximally entangled states |\E')mes = '^"^''^ IJ'"^); where 
d = {2j + 1) is the dimension of the Hilbert space. Utilizing = d'^/^, Eq. gives 
-Emax = log((i), which is the maximum value of the entanglement for the system [43| 144]. 
As an ansanz, the spin state {"^(t)) can be treated as a Guassian [T7] : 

. 1 (11) 

with its width (A J,) = (j/2)^/^^, determined by the number squeezing parameter ^^(t). 
Substituting Eq. ( ITTl) into Eq. (ITOl) . and replacing the discrete sum over m by an integral, 
we arrive at 

E{t)c:, hog [enjeit)], (12) 

which provides us analytical relation between the number squeezing and the two-mode 
entanglement. Considering = 1 for the initial CSS \j,j)x, we have E{0) = -Ecss — 
|log[e7rj]. Numerically, Hines et al. [13] have found that E'cssZ-^'max remains finite 
for large N. Our result shows -Ecss /-E'max — ^ | as — oo. In Fig. [2]^a) and (b), we 
plot time evolution of the entropy for finite A^ = 20 case. One can find that Eq. (fT2ll 
(red circles) agrees very well with exact numerical solution of Eq. fllOl) . In the initial 
stage, E{t) decreases from Eqss to its local minimum at tmm, due to the appearance 
of the MSS. In comparison with the CSS, the MSS approaches to localized Twin-Fock 
state [15]: |A^/2)i|A^/2)2 = |j,0), which exhibits g^^'^ = = E = 0. To confirm 
it, we plot evolution of probability distribution |cmP in Fig. (Hj^a). For the optimal 
coupling Q/x = 1.732 and A^ = 20, the system evolves into the MSS after a duration 
X^min = 0.284, which shows probability distribution peaked at m = 0. Quasi-probability 
distribution Q{6, 0) of the initial state is isotropic [Fig. [3]^c)], representing the minimal 
uncertainty relationship of a coherent state [7]. It becomes an elliptic shape due to the 
squeezing along and the anti-squeezing along Jy [Fig. [3] (d)]. 



I |2 1 
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Figure 2. (color online) The entropy of entanglement E (black solid lines) for = 20. 
The red circles are given by Eq. with numerical result of The initial value of 
E, Ecss = 5 log[e7r7V/2] = 3.21, and the maximal value £;,„ax = log(7V + 1) = 4.39. 
Other parameters arc the same with Fig. [1] 



Until now, we consider spin dynamics of Hamiltonian ([T]) with x > and the 
initial state \j,j)x- This state is a stable fixed point in phase space [12], so the phase 
coherence g^^\ the number squeezing and the two-mode entanglement E oscillate 
regularly with the same period. Using Eq. ([8]) and Eq. f|T2l) . it is possible to measure 
both and E by extracting g^^^ (i.e., the visibility) in atomic interference experiment 
[H [21 El m [5l [6l [31], [38] . There are two alternative experimental setups. One possibility 
is to trap the condensed ^^Na atoms in a symmetric double-well potential formed by 
atom chip [6], the other is the two-component BEG in an optical dipole trap [30]. In 
both cases, positive and large enough x ~ (c^ii + 0.22 — 2ai2) is required, which speeds 
up dynamics of the system such that the deleterious effects like atom losses can be 
neglected [25] . 

To proceed, let us consider another scenario: spin dynamics of Hamiltonian ([T]) 
with negative x case [16]. Now, the initial state \j,j)x corresponds to an unstable point 
at the separatrix [42], which leads to a more complex dynamics with a quite different 
characteristic time scale. From Eq. ([8]), we find that there is no number squeezing for 
negative x case [see also Fig. [Tj^c) and (d)]. Instead of preparing the MSS, the latter 
model can be used to generate Schrodinger's cat state likes, |\1/)cat = :^(b5 ~j) + IJ^J)) 
[47] . In Fig. [2](c) and (d), we plot time evolution of E{t) for = 20 case. Our results 
show that Eq. ( |T2i) follows exact results of E in the early stages of the evolution, then 
diverges as |\E') — >• |\1/)mes- This is because population distribution of |\I/)mes is no 
longer to be a Gaussian, and Eq. (fT2!) can not simulate the entropy. However, right-side 
of Eq. ([I2D, a monotonic function of is still useful to recognize the appearances of 
the cat state, which exhibits the largest number variance (AJ^)^^^^ = and ^'^^-^ = 2j. 
In real evolution, an approximate cat state with ~ 1.3j [Fig. [^d)] can be obtained 
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at a time Ixl^max ~ ln(8A^)/A^ for the optimal coupling = N/2 [17]. From right 

panel of Fig. [3l we also find that this state shows probability distribution |cmP peaked 
at m = ±j and maximal values of Q{6., 4>) pointed to the north and the south poles of 
the Bloch sphere. 




Figure 3. (color online) Probability distribution |c,„p (a)-(b), and quasi-probability 
distribution Qie,(t)) = | (6*, ^ (c)-(f), where the CSS |6',(/)) is given in Eq. Q. 

The red lines in (a) and (b) represents |cmP for the MSS at tn^n = 0.284x^^ and the 
cat state at tmax = 0.275|x|~^, respectively. The parameters in the left and the right 
panels are the same with Fig.[TJa) and Fig.[ljd). 

In summary, we have investigated spin dynamics of a symmetric BEC with repulsive 
interactions (x > 0) evolved from a coherent spin state \j,j)x- As main results of our 
work, we find analytical expressions of the number-squeezing parameter ^"^ and the 
entropy of entanglement E as, Eq. ([8]) and Eq. ( fT2l) . Both of them can be, in principle, 
at least, measured by extracting the phase coherence g^^\t) (i.e., the visibility) in atomic 
interference experiments. For the case of attractive interactions (x < 0), thought there 
exists no number squeezing, the squeezing parameter or Eq. f|T2|) is still useful to 
recognize the appearance of Schrodinger's cat state. 
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